Redshift Evolution of the Low-Column Density Lya Forest 



Esther M. Hu 1 

1 Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 

66 ' 

0^ ■ Abstract. Knowledge of the evolution of the low column density end of the Lya 

G\ | forest (12.8 < log N m < 16.0 cm 2 ) can provide substantial insight for models of 

structure formation and studies of the evolution of the ionizing background. Here 
we describe the evolution of the b values and column densities of forest clouds over 
the redshift range 2 < z < 4. While the distribution function for clouds with column 
densities iVm < 3 x 10 14 cm" 2 can be well described by a single invariant power 
law ht with slope ~ — 1.5 over this range, at higher column densities a break in this 
distribution function is seen, with a deficiency of high column density clouds relative 
^ ' to the single power-law fit, and with the break point moving to lower column densities 

for clouds at lower redshifts. The median and cut-off (minimum) b values are seen 
■ to decrease at higher redshifts. 
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1 Introduction 



Recent CDM-based numerical simulations p5|,||4j],|lq] have shown that the 
Lya forest clouds are a natural consequence of structure formation in the in- 
tergalactic gas, and a general picture has emerged within which their 

properties may be understood. In particular, the low column density (12.8 < 
log Njji < 16.0 cm -2 ) forest corresponds to underdense or slightly overdense 
^ | regions in the IGM, which are easily modeled or analytically simulated, and 

which contain a substantial fraction of the baryonic material in the intergalac- 
tic gas. Studies of this component over as wide a redshift range as possible 
■ can provide information which: (1) constrains models for structure formation, 

(2) provides normalization for studies of Hell, and (3) addresses the evolution 
of the shape and normalization of the ionizing background radiation. 

The advent of large, homogeneous data sets from the Keck 1 10-m telescope 
obtained with the HIRES spectrograph at high resolution (R ~ 36, 000) and 
high signal-to-noise || , || has made it possible to explore the evolution of the 
low column density forest in some detail. The present analysis uses data from 
the spectra of five quasars (Q1623+268, Q1700+643, Q0014+813, Q0302-003, 
Q1422+231) which span a redshift range from z em = 2.521 — > 3.620, and 
whose forest clouds have been analyzed over a range from z a bs = 2.17 — > 3.51. 
This wavelength region was chosen to sample the 'quiet' forest, and to avoid 
the region of the proximity effect, and in the case of two partial Lyman limit 
systems, to avoid the region within ±1000 km s" 1 of these systems. We have 
extended the redshift range treated in this work to z ~ 4 using the results 
(3.40 < z abs < 4.00) of Lu et al. 1996 [pi, who employ a different fitting 



algorithm but use similar modeling analysis to correct for selection effects 
and incompleteness. Details of the profile fitting can be significant |0|, as 
is discussed in more detail by Kim et al. JlO| in these proceedings, and a 
consistent analysis over the full redshift range using a single procedure should 
ultimately be carried out, but the Lu et al. data seems fully consistent with 
the extrapolated trends from z = 2.2 — > 3.5. 

2 The Differential Density Distribution Function 

The spectrum of HI column densities can be characterized by a power law 
with slope ~ —1.5 over a wide range of column densities (Figure I) but with 
a marked deficiency of clouds relative to this fit in the region 10 143 < Nbi < 

10 16.3 cm -2 
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Figure 1: Spectrum of HI column densities at z ~ 3. The best power-law 
fit (oc iV^j 1 ' 46 ) to the number of clouds per unit column density obtained for 
the weak forest clouds (filled squares) at (z) = 2.85 by Hu et al. 1995 [|| is 
shown extrapolated through the points of Petitjean et al. 1993 [[lJJ (diamonds, 
with la error bars) out to column densities of I0 22 cm~ 2 using their notation 
for this function, and provides a surprisingly good description over nearly 10 
orders of magnitude in column density. The deficiency of clouds noted above 
Nm > 10 143 cm -2 is a real effect which cannot be attributed to mis-estimates 
in the column density due to saturation effects. 




Figure 2: (Left Panel): The distribution of HI column densities by redshift 
bin, with ler error bars. The solid line shows the power-law fit (<x N-^f ) 
shown in Figure 1 and the dotted line shows the single power-law fit obtained 
by Lu et al. for the (z) = 3.7 data (open circles). At Nm < 10 14 cm -2 , 
the distribution functions at different redshifts look very similar, however, at 
^Vhi > 10 14 cm -2 , the deficiency in the number of forest clouds compared 
to the expected power-law distribution depends on redshift. The point of 
deviation from the power law moves to smaller Nm as redshift decreases, 
indicating the rapid evolution in the higher Nm forest clouds. (Right Panel): 
Incompleteness corrected differential density distribution functions compared 
for (z) — 2.31 and (z) = 3.35 systems. The reference lines show the power-law 
fits of the left-hand panel. 



Kim et al.'s 1997 || analysis of the differential density distribution function 
(DDF) by redshift interval clearly shows an evolution in the departure from 
the power-law fit (Figure 2). In the left-hand panel of Figure 2 the (z) = 3.70 
data jjj] show only a slight steepening (dotted line) at the higher column 
densities, but by (z) — 3.35 a deficiency in the number density of clouds at 
Nm > lO 14 8 cm -2 can be seen, while clouds at (z) = 2.85 and (z) — 2.31 
show deviations from the power-law extrapolation at A/hi > 10 143 cm~ 2 . This 
departure from a power law can be easily seen in the right-hand panel of 
Figure 2, where the distribution functions of the (z) — 2.31 (diamonds) and 
(z) = 3.35 (triangles) absorbers, now corrected for incompleteness, are shown. 
There are no clouds with TVhi > 10 15 2 cm~ 2 at (z) = 2.31 despite the fact that 
the redshift path is 1.5 times larger than at (z) = 3.35, so the deficiency of high 
column density clouds in the lower redshift range is not due to observational 
selection. We conclude that we are seeing a rapid evolution with redshift in 
the number density of absorbers at high column density, with the break in 
the density distribution function strengthening and migrating to lower column 
densities as redshift decreases. 




Figure 3: N(z) vs. redshift for the present data (filled circles) and the HST ob- 
servations (open squares] open diamonds). The solid line shows the maximum 
likelihood fit to the HST data and the high column density Keck data (slope 
7 = 2.15). The dotted line shows the fit to the Keck data only, with a slope of 
7 = 2.78. Triangles show the forest clouds over just the column density range 
N m = 10 131 - 10 14 cm" 2 , with a long dashed line fit with 7 = 1.19. The la 
error bars and redshift binning are indicated for each point. 



Turning to the number density evolution of Lya clouds with redshift, where 
the number density of clouds per unit redshift per line of sight is expressed 
as: N(z) = Nq(1 + z) 7 , we can compare the Keck data on the high redshift 
systems with the low redshift systems studied with HST , || , while keeping 
in mind that N(z) will be sensitive to the adopted threshold Nm (a reflection 
of the departure from the single power-law description at Nm > 10 143 cm -2 ). 
Figure 3 shows a plot of log N (z) vs. log (1 + z) giving the maximum likelihood 
fit to the HST data and the high column density Keck data (7 = 2.15). The 
adopted column density range used for the fit, Nm = 10 1377 — 10 16 cm~ 2 , 
corresponds approximately to the conventional equivalent width threshold of 
W = 0.32 A for comparisons with the lower redshift data ]lj,|^],||. By con- 
trast, the best fit to the slope of N(z) for the low column density systems in 
the Keck data is 7 = 1.19. The definition of N(z) yields 7 — 1 for q = 
and 7 = 0.5 for q = 0.5 in the case of non-evolving clouds, and the slopes 
for the data in different column density ranges in Figure 3 can be understood 
in terms of the evolution of the break in the DDF: the number density of the 
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Figure 4: Observed (dotted histogram) and simulated (solid histogram) b dis- 
tributions. For the (z) — 3.35 clouds, the low cut-off b value is ~ 17 km s _1 . 
For the (z) = 2.31 clouds, there is no single best-fit assumed distribution 
for b. Assuming that most narrow lines are unidentified metal lines towards 
Q 1700+643 and ignoring the lines with b < 22 km s _1 , the cut-off b of 
24 km s _1 which is shown in the figure best fits the data. 



higher column density systems (7 = 2.15) is evolving rapidly, while the num- 
ber density of the lower column density systems (10 131 < JVjji < 10 14 cm~ 2 ; 
7 = 1.19) is not. 



3 The Evolution of b- Values 

A second interesting issue is the evolution of line-widths, characterized by the 
Doppler parameter b obtained from profile fitting, whose behavior reflects the 
thermal/ionization history of the IGM ||. Hu et al. 1995 |J found a typical 
mean b of 28 km s _1 with a minimal b c = 20 km s _1 from the analysis of 
(z) = 2.85 systems. At (z) = 3.70 Lu et al. (l2) report a mean b of 23 km s _1 
and cutoff b c = 15 km s , showing a trend to narrower absorbing components 
at higher redshift. 

In order to quantify these descriptions, it is important to run model sim- 
ulations to understand the nature of blending, selection, and incompleteness 
effects on the analysis of the observed data. The procedure used ||,[|l^] is 
to input a truncated Gaussian distribution for the Doppler parameter values, 
generate a series of simulated artifical spectra, and then apply the standard 
profile-fitting methods to match the recovered column densities and b values 
to the input spectra and to the observations. 



Both the cutoff and median b values appear to increase with decreasing 
redshift — the cutoff b value's increase from 15 km s _1 at (z) = 3.70, to 
17 km s _1 at (z) = 3.35, to 22 km s" 1 at (z) = 2.85, to 24 km s _1 at (z) = 2.3 
presumably reflecting an increase in the turbulent and thermal broadening 
over this period. 



4 Summary 

In broad outline the cloud evolution described here probably follows naturally 
from the decreasing IGM density with decreasing redshift, which drives down 
the column density at which shocking occurs ]13|| , and the increasing level of 
structure which results in more turbulent components. However, reproducing 
the results in detail will require full comparison with a wide suite of numerical 
models, which may provide useful constraints on the history of the ionization 
and heating of the IGM. 
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